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Motivation: 

Large-Eddy  Simulation  (LES)  Challenges 
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Algorithm  Comparisons 

-  identical  subgrid  modeling 

-  differences  reside  in  numerics 
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Ref:  2013  -  Cocks,  Sankaran,  Soteriou,  “Is  LES  of  reacting  flow  predictive?  Part  1:  Impact  of  Numerics” 


Need  to  determine  BEST  discretization  scheme  for  Reacting  LES  | 
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Approach 


-  Investigate  dissipation  and  dispersion  characteristics  of 
schemes 

-  tied  to  solution  accuracy 

-  use  Von  Neumann  Stability  Analysis 
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-  Schemes  to  investigate: 
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-  Standard  Collocated  Grid  ^ 

-  Standard  Staggered  Grid 
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FIG.  2.  Staa-eered  arid  system. 
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Von  Neumann  Analysis 


— +  — =  0  IDEulerEqns  — >  *^=-+  A*^=-  =  0  with^4  =  dE / dQ 

dt  dx  dt  dx 


Eigenvalues  of  the  amplification  matrix 
specify  growth  factor  and  phase  errors. 
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Staggered  Grid  Scheme/  Quasi-Linear  Form 
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Stability  Analysis 
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Stability  Analysis 
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Test  Cases 


ID  Duct 

-  Periodic  BC’s 
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Saw-tooth  i.c.  Sinusoidal  i.c. 
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Particle  Wave  High  Frequency  Behavior 


Runge  Kutta  Scheme! 
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Growth  Factor 
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Acoustic  Wave  High  Frequency  Behavior 


Runge  Kutta 


Mach  =  0.28 

_ distrub  =  0.01% 

Staggered: 

-  DAMPED  with  mild 
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Solution  @  i: 


phase  error 


Collocated: 

-  stationary  &  un 
damped 
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Effect  of  Boundary  Conditions 


A2At 

CFL  =  1  = 

Ax 

Mach  =  0.28 
distrub  =  0.01% 


Pressure  at  10T3 
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characteristic  variables  at  10T3 


-  high  frequency  in  Collocated  pressure  solution 
-  lack  of  acoustic  damping 
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Kinetic  Energy  Preservation  (KEP) 


“in  computations  of  turbulent  flow  fields,  dissipative  errors  show  up  at  the  level  of  kinetic  energy” 
(Mahesh  2004) 

Robust  at  inviscid  limit  (Re  ■>  oo) 


Incompressible  Flow: 
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Compressible  Flow: 


K  =  y2  Uj2  bounded  and  constant  at  inviscid  limit 
KEP  schemes  satisfy  secondary  equation  discretely 
Richtmeyer  &  Morton  (1967) 

Arakawa  (1966) 
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Discrete  analogue  seeks: 

Accurate  transport  of  KE  -^accurate  physical  transfer  of  energy:  E  =  KE  +  Uint 
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KEP:  Applied  to  ID  Euler 


(Collocated  Grid)! 


compare  Crank-Nicolson  (CN)  with  Fully  KEP  scheme  (F-KEP) 
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discrete  secondary  equation  satisfied  to  machine  zero  if  KEP 
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Evaluating  KEP:  u2 
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Crank-Nicolson  Scheme 
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F-KEP  always  secondary  conservative 
CN  KEP  for  low  compressibility  effects 
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Going  Forward 


•  Key  Questions: 

•  What  is  the  advantage  of  kinetic  energy  preservation  for  LES? 

•  Does  it  help  minimize  or  eliminate  the  need  for  artificial  dissipation? 

•  What  about  the  relative  importance  of  dispersion  errors? 

•  Implement  Merkle’s  generalized  KEP  schemes 

•  Formulated  for  both  staggered  and  collocated  schemes 

•  Major  advantage  is  that  it  is  KE  preserving  for  the  scalar  transport  as  well 

•  Can  we  minimize  or  eliminate  the  need  for  artificial  dissipation  for  scalars? 

•  Extend  schemes  to  multi-dimensional  code 

•  Apply  to  non-reacting  and  reacting  LES  computations 
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